Extended Phonon Collapse and the Origin of the Charge-Density- Wave in NbSe2 
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We report inelastic x-ray scattering measurements of the temperature dependence of phonon 
dispersion in the prototypical charge-density- wave (CDW) compound NbSe2. Surprisingly, acoustic 
phonons soften to zero frequency and become overdamped over an extended region around the 
CDW wavevector. This extended phonon collapse is dramatically different from the sharp cusp in 
the phonon dispersion expected from Fermi surface nesting. Instead, our experiments combined with 
ab initio calculations, show that it is the wavevector dependence of the electron-phonon coupling 
that drives the CDW formation in NbSe2 and determines its periodicity. This mechanism explains 
the so far enigmatic behavior of CDW in NbSe2 and may provide a new approach to other strongly 
correlated systems where electron-phonon coupling is important. 
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The origin of CDW order is a long-standing problem 
relevant to a number of important issues in condensed 
matter physics, such as the role of stripes in cuprate 
superconductivity [ij and charge fluctuations in the colos- 
sal magnetoresistive manganites 2] . Static CDW order, 
i.e., a periodic modulation of the electronic density, re- 
flects an enhancement of the dielectric response of the 
conduction electrons at the CDW wavevector, (\cdw, 
but it has long been known that it is only stabilized by 
a coupling to the crystal lattice^, Transitions into 
the CDW phase on lowering the temperature are accom- 
panied by a softening of an acoustic phonon at (\cdw 
to zero frequency at Tcdw where it freezes into a static 
distortion 5] and evolves into the new periodic (often in- 
commensurate) superstructure. Chan and Heine derived 
the criterion for a stable CDW phase with a modulation 
wavevector q asjJl 
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where rjq is the electron-phonon coupling associated with 
a mode at a frequency of ujq, Xg is the dielectric response 
of the conduction electrons, and Uq and Vq are their 
Coulomb and exchange interactions. Although both sides 
of this inequality are essential in stabilizing the CDW 
order, the common assumption is that the modulation 
wavevector, qcDW, is determined by the right-hand side, 
i.e., by a singularity in the electronic dielectric function 
Xq at q_cDW = 2A:f, where kp is the Fermi wavevector. 
In the case of NbSe2 , it was proposed that such a singu- 
larity resulted either from direct Fermi surface nesting at 
Qc DW 0, @l from the presence of saddle points near 
the Fermi surface connected by qcniy [ll- However, this 
has been challenged by a recent density-functional-theory 



(DFT) calculation that correctly predicted a CDW in- 
stability, but without singularities in Xq S • Further- 
more, several angle-resolved photoemission spectroscopy 
(ARPES) experiments have found no clear evidence of an 
energy gap at the nesting wavevectors opening at Tcdw 
jlQi-il2i] . As a consequence, the elegant and intellectually 
compelling picture of CDW formation driven by Fermi 
surface nesting has been called into question 's^. The al- 
ternative possibility is that the CDW transition in NbSe2 
may instead be driven by the left-hand side of equation 
1, i.e., the wavevector dependence of the electron-phonon 
coupling r]q. A direct test of this conjecture by phonon 
spectroscopy is the subject of this letter. 

We have investigated discrepancies with the Fermi 
surface-nesting scenario by measuring the energies and 
linewidths of phonon excitations in NbSe2 as a function 
of temperature, comparing both to ab initio calculations. 
In particular the phonon linewidth, which is proportional 
to the product of the electron-phonon coupling rjq and the 
electronic joint density-of-states (JDOS), which in turn 
is related to Xg: provides detailed information related to 
the Fermi surface. Our results demonstrate that the soft 
phonon physics in NbSe2 does indeed deviate strongly 
from the conventional picture based on nesting and 
that the electron-phonon coupling is primarily responsi- 
ble for determining the wave vector of CDW order. 

Earlier inelastic neutron scattering investigations of 
the phonon dispersion in NbSe2 were limited by weak 
scattering intensity due to the small sample size and by 
relatively poor wavevector resolution 13l |14| . To over- 
come these constraints, we utilized high-resolution in- 
elastic X-ray scattering (IXS), which allowed us to ob- 
tain measurements of the entire dispersion of the soft 
mode branch in a wide temperature range. We used a 
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FIG. 1. (color online) Temperature dependence of the soft 
phonon mode and the charge-density-wave superlattice peak 
near qcDVK ~ (0.329,0,0). (a)-(d) Energy scans at q = 
(3 - ft, 0,0), h = 0.325, for temperatures 8K < T < 90K. 
Solid (red) lines are fits consisting of damped harmonic oscil- 
lators (inelastic) and a pseudo-voigt function (elastic) (blue 
dashed lines), (e) Intensity of the charge-density- wave super- 
lattice peak for T < 120 K. The inset shows the phonon and 
superlattice peak intensities just above Tc. (f)(g) Phonon fre- 
quency u)q and critical damping ratio F /ujo of the soft phonon 
mode, respectively, at q = {h, 0, 0) with h = 0.3 (circle), 0.325 
(triangle) and 0.35 (square). The solid line in (f) is a power 
law fit of the form ((T - Tc)/Tcf yielding 5 = 0.48 ± 0.02. 



high-quality single crystal sample of about 50 mg (2 x 
2 X 0.05 mm^) with a Tcdw of 33 K determined from 
the temperature dependence of the superlattice reflec- 
tions (Figure [J) in agreement with previous results [13[. 
All IXS experiments were carried out on the XOR 30- 
ID HERIX beamline at the Advanced Photon Source, 
Argonne National Laboratory. Data were fitted with 
damped harmonic oscillator functions convoluted with 
the experimental resolution. For more details on the in- 
strumental setup and data analysis see the supplemen- 
tary material. Here, we focus on the longitudinal acoustic 
phonon branch dispersing in the crystallographic (100) 
direction and crossing the CDW wave vector q_cDW — 
(0.329, 0, 0)[l3|. Figure [T] shows the temperature depen- 
dence of a soft phonon mode at q^ki = (0.325,0,0), 
close to the CDW wavevector, qcDW = (0.329,0,0). At 
T = 90K (Fig. [T]i), the soft phonon at an energy of 
LUq = 4.5 meV has nearly equal intensity to the second 
phonon branch at 10 meV. Upon cooling, the intensity 




-10 -5 5 10 15 
energy (meV) 



10 -5 5 10 15 
energy (meV) 



FIG. 2. (color online) Wavevector dependence of the soft- 
phonon at T = 33 K. Energy scans at Q = (3 — h,0,l), 
h = 0.275 — 0.375. Solid (red) lines represent the total fit 
result consisting of a damped harmonic oscillator functions 
(inelastic) and a pseudo-voigt function (elastic) (blue dashed 
lines) . 



of the upper branch is suppressed due to the Bose factor, 
whereas the intensity of the soft phonon is enhanced by a 
factor 1/iOq in the cross section as its frequency, Wg, is re- 
duced. At T = Tcdw, the elastic superstructure peak of 
the CDW phase dominates the spectrum (Fig. [Hd), but 
we can still distinguish the critically damped phonon as a 
broad peak beneath the narrow elastic CDW peak. Well 
inside the CDW phase, the elastic superlattice refiection 
was too strong for any inelastic scattering to be observed 
at cicDW and these low temperature data (Fig. [l}i) were 
used to determine the shape of the elastic scattering. 

Figure [T]3 shows that the integrated intensity of the 
CDW superlattice peak measured at /i = 0.325, which 
is within the momentum- resolution of qcDVK, increases 
rapidly below Tcdw = 33 K, in good agreement with 
previous neutron diffraction data taken on crystals from 
the same growth batchfis']. Above Tcdw, the elastic in- 
tensity due to diffuse scattering from the sample is very 
small, which implies that our sample had very little struc- 
tural disorder. It stays low until very close to Tcdw 
(see inset to Fig. [T^), where a weak elastic "central" 
peak consistent with low energy critical fluctuations ap- 
pears. The phonon frequency at qcDW softens on cooling 
(Fig. Hf) following a power law a;,(T) = ((T - T^)/T^Y 
with S = 0.48 ± 0.02, the value predicted by mean- 
field theory Q. As the phonon softens, the damping in- 
creases and the phonon becomes critically damped, i.e. 
F/wq — 1, at Tcdw (Fig. HT). Remarkably, we observe 
the same power law behavior not only at ciCDW, but also 
at h — 0.3 and 0.35 which are far outside the experimen- 
tal resolution from qcDW and where the elastic peak is an 
order of magnitude weaker relative to the phonon inten- 
sity. Moreover, the phonon frequencies at these wavevec- 
tors also become indistinguishable from zero at Tcdw 
(Fig. [33,c) the same as at cicdw- This means that 
the phonons are critically damped over a large range 
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FIG. 3. (color online) Experimentally obtained dispersion 
and damping ratio of the soft phonon branch in NbSe2 at four 
temperatures 8K < T < 250 K. Plotted are (a) the frequency 



of the damped harmonic oscillator (DHO) 



and (b) the damping ratio T/ljo- Lines are guides to the eye. 
Note that phonons at /i = 0.325, 0.35 and T = 8K were not 
detectable due to strong elastic intensities. The inset in (b) 
shows the experimentally observed damping F of the damped 
harmonic oscillator (symbols) and scaled DFPT calculations 
(see Fig. g]) of 27 (lines, offset 0.7meV) with a = 0.1 eV 
(black) and 1 eV (red). 
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FIG. 4. (color online) Ab initio calculation for the soft phonon 
mode along the crystallographic (100) direction in NbSe2. 
Shown are the calculated (a) dispersion, (b) 27 (the con- 
tribution of the electron-phonon interaction to the phonon 
linewidth (FWHM)), and (c) the electronic joined density 
of states (JDOS). Calculations were done with five different 
Gaussian broadenings 0.1 eV < cr < 1.0 eV (see text). Note 
that calculated imaginary frequencies are shown as negative 
roots of the square phonon frequencies. Lines are guides to 
the eye. 



of momentum transfer from h = 0.3 to h ~ 0.35. Go- 
ing further away from cicDW with the same step size, 
Ah = 0.025 r.l.u., the soft phonon branch is well sep- 
arated from zero energy (Fig. [2k, d). Figure [3] shows 
the full dispersion and damping ratio of the soft-mode 
phonon branch. A broad dispersion anomaly is already 
evident at 250 K in agreement with previous neutron 
scattering measurements performed only at 300 K^. 
This anomaly deepens considerably upon cooling to 50 K, 
where we also observe a strong increase in the damping. 
Finally, upon cooling to TcdW: the frequencies reach zero 
and the phonons become critically damped over an ex- 
tended range of wavevectors. At T = 8K, well below 
TcDWi we find hardened frequencies and reduced damp- 
ing, similar to the ones observed at T = 50 K. However, 
the soft mode was not resolvable at h — 0.325 and 0.35 
due to strong elastic scattering from the CDW peak (Fig. 
[Iji) . At these temperatures the Bose and 1 /uj factors sup- 
press the phonon intensity and the measurements become 
increasingly difficult. 

The q-dependence of the phonon softening shown in 
Figure[3]is in marked contrast to the sharp, cusp-like dips 
that normally characterize Kohn anomalies at 2\ip due 
to Fermi surface nesting^, [T^. For example in the CDW 
compound ZrTes, where the importance of Fermi surface 
nesting has been clearly established [I^, phonon frequen- 
cies are only renormalized over a q-range of 0.085 A^^ 
and soften to zero at a single wavevector using a step size 
of 0.04A~^[i|. In contrast, in NbSe2, we find that the 
phonon renormalization extends over 0.36A-\ or over 



half the Brillouin zone, and the critically damped region 
extends over 0.09 A^^, about twice the experimental res- 
olution. This behaviour clearly rules out a singularity 
in the electronic response in NbSe2 and suggests that 
the CDW is determined by the wavevector dependence 
of the electron-phonon coupling riq, as proposed by Jo- 
hannes et al.Q. A broadened or even flat-topped sus- 
ceptibility due to imperfect nesting caused, e.g., by the 
c-axis dispersion of the electron bands, could also lead to 
a renormalization of the phonon dispersion over a larger 
range of wavevectors, but it is unlikely that it spans over 
half of the Brillouin zone. 

In order to elucidate the microscopic mechanism be- 
hind the CDW phase transition in NbSe2, we compare 
our experimental results to detailed phonon calculations 
based on density functional perturbation theory (DFPT) 
performed with the crystal structure at T > Tcdw- This 
is a zero temperature technique, in which structural in- 
stabilities show up as imaginary phonon frequencies (see 
supplementary material). Because of the finite momen- 
tum mesh used in the DFPT calculations a numerical 
smearing, a, of the electronic bands is necessary to com- 
pare the calculations with experiment. The effect of a is 
analogous to a thermal smearing of the electronic struc- 
ture, so it has been used in previous work to simulate the 
effect of temperature 15|, Il7| . However, a naive interpre- 



tation, i.e., a = 2.12 * A;BT[l5j, results in temperatures 
at least one order of magnitude too large with respect to 
experimental observations 17j • We believe that the cal- 
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culated temperatures disagree quantitatively with exper- 
iment because the calculations ignore the thermal disor- 
der of the lattice, which can have profound effects on the 
phonon excitations . Elucidating these issues is an im- 
portant direction for future research beyond the scope of 
the present work. We merely note that for NbSe2 a com- 
parison between theory and experiment indicates that 
values of 0.1 eV < cr < 1 eV produce results that are con- 
sistent with a temperature range of 30K < T < 300 K. 

Figure |4] summarizes the calculations, showing the 
calculated soft-phonon dispersion, linewidth, and elec- 
tronic joint density of states (JDOS). Imaginary phonon- 
frequencies are represented in Fig. through the nega- 
tive roots of the absolute value, e.g., as 'negative' phonon 
frequencies. These occur in the calculated longitudi- 
nal acoustic phonon branch for a > 0.18 eV over an 
extended range of wavevectors (Fig. |4^) in agreement 
with previous studies [3, and in a qualitative agreement 
with the observed breakdown of the phonon dispersion. 
Similarly, the contribution to the phonon linewidth from 
the electron-phonon interaction, 27 (Fig. |41d), shows a 
strong enhancement over the same extended range of 
wavevectors. 27 is proportional to the product of the 
electron-phonon coupling and the electronic joint den- 
sity of states (JDOS). Since the latter shows negligible 
wavevector dependence (Fig. |4j:), the enhancement of 
the phonon linewidth observed in both experiment and 
theory is entirely due to a strong wavevector dependence 
of the electron-phonon coupling. Moreover, this range 
of strongly enhanced electron-phonon coupling is iden- 
tical to the range over which the phonon softens. In 
contrast, the real part of the susceptibility only shows a 
much broader and only very shallow peak 8]. This leads 
us to conclude that the observed q-dependence of the 
phonon self-energy is entirely due to the electron-phonon 
coupling and that the CDW wavevector in NbSe2 is in- 
deed determined by the left-hand side of equation [TJ 

Since our measurements demonstrate that Fermi sur- 
face nesting does not play a role in the CDW phase tran- 
sition in NbSe2, the electronic states serve only to pro- 
vide an elevated dielectric response, with the modulation 
wavevector entirely determined by the coupling between 
electronic and vibrational dynamics. Previous studies of 
chromium[19] and ruthenium [20| have shown that ma- 
trix elements can indeed be sharply wavevector depen- 
dent and also produce dips in the phonon dispersions, 
although in these compounds phonons do not soften to 
zero energy. Our work provides direct evidence that the 
same effect can drive the structural instability in a CDW 
compound. This result naturally explains why electronic 
probes do not find strong nesting or a clear gap at q^cDW 
in NbSe2. Our results have implications for many other 
strongly correlated systems. In particular, CDW correla- 
tions in the form of stripes and/or checkerboard patterns 
have been linked to the emergence of unusual states and 
physical properties, such as colossal magnetoresistance 



in the manganites[21] and the pseudogap state in the 
cuprates[2J]. Indeed, the observation of phonon anoma- 
lies in manganites at the wavevector of the checkerboard- 
type order [2 3j and anomalies observed in La2-2:Sr2:Cu04 
at the stripe ordering wavevector j24j demonstrate that 
strong electron-phonon coupling could be important in 
these materials as well. 

In conclusion, we reported inelastic x-ray measure- 
ments of the temperature dependence of a longitudinal 
acoustic phonon in NbSe2 involving the CDW soft mode. 
We observe an extended region in q with overdamped 
phonons at the CDW transition temperature. A de- 
tailed comparison to lattice dynamical calculations via 
DFPT shows that in NbSe2, the periodicity, i.e. qcDW, 
of the CDW ordered state is determined entirely by the 
wavevector dependence of the electron-phonon coupling. 
This is in stark contrast to the standard view that a di- 
vergent electronic response defines c^cdw ^-nd is evidence 
that a wavevector dependent electron-phonon coupling 
can drive a structural phase transition. 
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